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Abstract 

In previous papers and letters on quantum amplitudes in black-hole evaporation, a boundary-! 
value approach was developed for calculating (for example) the quantum amplitude to 
have a prescribed slightly non-spherical configuration of a massless scalar field <f> on a 
final hypersurface Sf at a very late time T , given initial almost-stationary spherically- 
symmetric gravitational and scalar data on a space-like hypersurface Ej at time t = . 
For definiteness, we assumed that the gravitational data are also spherically symmetric on 
Tip. Such boundary data can correspond to a classical solution for the Einstein/scalar 
system, describing gravitational collapse from an early low-density configuration to a 
nearly- Schwarzschild black hole. This approach provides the quantum amplitude (not 
just the probability) for a transition from an initial to a final state. For a real Lorentzian 
time-interval T , the classical boundary- value problem refers to a set of hyperbolic equa- 
tions {modulo gauge), and is badly posed. Instead, the boundary- value approach of the 
previous letters and papers requires (following Feynman) a rotation into the complex: 
T — > |T|exp(— id), for < 9 < tv/2, of the time-separation-at-infinity T. The classical 
boundary-value problem, for a complex solution of the coupled nonlinear classical field 
equations, is expected to be well-posed for < 9 < %/2 . For a locally-supersymmetric 
Lagrangian, containing supergravity coupled to supermatter, the classical Lorentzian ac- 
tion S'dass, a functional of the boundary data (which include the complexified T), yields 
a quantum amplitude proportional to exp(i5 , c i aS s), apart from possible loop corrections 
which are negligible for boundary data with frequencies below the Planck scale. Finally 
(still following Feynman), one computes the Lorentzian quantum amplitude by taking the 
limit of exp^Sdass) as 9 — > 0+ . In the present paper, a connection is made between 
the above boundary- value approach and the original approach to quantum evaporation in 
gravitational collapse to a black hole, via Bogoliubov coefficients. This connection is devel- 
oped through consideration of the radial equation obeyed by the (adiabatic) non-spherical 
classical perturbations. When one studies the probability distribution for configurations 
of the final scalar field, based on our quantum amplitudes above, one finds that this dis- 
tribution can also be interpreted in terms of the Wigner quasi-probability distribution for 
a harmonic oscillator. 

1. Introduction 

In previous letters and papers [1-3] , we calculated the quantum amplitude for a given 
configuration of the massless scalar field </> on a space-like hypersurface Tip at a very late 
time T, measured at spatial infinity, given (for simplicity) that both the gravitational and 
scalar field configurations (hij, 4>)i on an initial hypersurface Ej are exactly spherically 
symmetric, and that the gravitational data hijF on Ej? are also exactly spherically symmet- 
ric. Here, hij = g%j denotes the intrinsic spatial metric = 1,2,3), while <7 M „ denotes 
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the 4-dimensional space-time metric (fj,, v = 0,1,2,3). Mathematically, this calculation 
involves analytic continuation of the quantum amplitude into the lower complex T-plane. 
In [2,3] on the 'complex approach' and the 'spin-0 amplitude' in black-hole evaporation, 
this was phrased in terms of a rotation of T into the complex: T — > |T| exp(— i9), where 
< 9 < 7r/2. Provided that the appropriate system of coupled Einstein/massless-scalar 
classical field equations is strongly elliptic [4], up to gauge, for < 9 < tt/2 , the complex 
classical boundary- value problem, corresponding to the above quantum amplitude, should 
be well posed. In the extreme case 9 = n/2 , the boundary data have been rotated so as to 
provide a boundary-value problem for a real Riemannian (positive-definite) 4-metric , 
coupled to a real scalar field (p. The field equations are 'elliptic modulo gauge', for which 
one expects good behaviour of the boundary-value problem, including analyticity of the 
solutions. 

It is assumed throughout that the original bosonic theory, such as Einstein gravity 
with a minimally-coupled massless scalar field, is embedded in a locally-supersymmetric 
theory containing supergravity with supermatter. In this case, the smallest such theory 
[5] contains N = 1 supergravity with its spin-3/2 gravitino, a complex scalar field and its 
massless spin-1/2 partner. Only for locally-supersymmetric theories does one expect that 
the loop terms Aq, Ai, A 2 , . . . in a semi-classical expansion of the quantum amplitude 
Amp ~ (Aq + HAi + h 2 A 2 + . . . ) exp(— Ib/K) will be finite [6-8]. Here, for our bosonic 
boundary data above, Ib is the classical 'Euclidean action' of the boundary-value solution 
of the coupled Einstein and bosonic-matter field equations. For simplicity, in our example, 
at present the only non-zero boundary data allowed are a real Riemannian 3-metric hij 
and a real scalar field (p . Note that, for dimensional reasons, the loop corrections will only 
be significant for boundary data which contain frequencies at and beyond the Planck scale. 

In [3] , the quantum amplitude was then evaluated semi-classically for 9 > , and 
finally the limit of the amplitude was taken as 9 — > 0+ , to obtain the Lorentzian amplitude 
(that is, the amplitude for T real), following Feynman's +ie prescription [2]. In particular 
(Eq.(4.16) of [3]), we computed the imaginary part of the (necessarily) complex spin-0 
classical action by taking the contribution from the poles along the real-frequency axis. The 
discrete set of real frequencies a n = mr/\T\ , for n = 1, 2, 3, . . . , was found to dominate 
in late-time field configurations. So far, however, we have made no explicit mathematical 
connection with the usual theory of black-hole evaporation [9-12]. We rectify this omission 
in this paper, by considering Bogoliubov transformations between (real) massless spin-0- 
field modes propagating in a slowly-varying Schwarzschild-like solution, such as appears 
in Sec. 2 of [3] for adiabatic modes. Such Vaidya-like approximate classical solutions will 
subsequently be treated in more detail [13]. 

The original derivation of black-hole radiance involved the calculation of Bogoliubov 
transformations between initial and final quantum states. Because of the time-dependence 
of the gravitational collapse, an initial basis of (linearised) positive- and negative- frequency 
massless-scalar modes (for example) differs from a final basis. By completeness, there exists 
a (Bogoliubov) transformation between the initial and final bases. When this transforma- 
tion mixes positive and negative frequencies, as in the black-hole case, one interprets the 
result physically in terms of particle creation. 

The connection between the approach used in [3] and the approach using Bogoliubov 
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coefficients arises from the radial equation satisfied by the (adiabatic) perturbations. The 
real (spin-dependent) potential in the radial equation vanishes sufficiently rapidly at spa- 
tial infinity that the radial functions near infinity (for each spin s) are superpositions of 
complex exponentials with complex coefficients; see Eq.(3.3) of [3] for the case s = . In 
the approach of [3] , regularity conditions for the fields at small radius imply that the radial 
functions are real. An early-time basis of modes, corresponding to waves incoming from 
past null infinity, can be related, by completeness, to a late-time basis, which is a superpo- 
sition of waves outgoing at future null infinity. It is then possible to relate, on a common 
late-time space-like hypersurface, the Bogoliubov coefficients to the complex coefficients 
above, which arose in the boundary-value method, since the Bogoliubov coefficients are 
independent of space-like surface. The usual (Bogoliubov) approach gives a density ma- 
trix and probabilities for final configurations, but not the phase information present in a 
quantum amplitude. 

Working instead with the boundary- value approach [3] , we have, as seen above, a pure 
state, and can evaluate probabilities as |amplitude| . Thus, our approach exhibits a large 
qualitative difference, as compared to the usual approach. Correspondingly, there should 
be some differences between the physical predictions based on the Bogoliubov approach 
and on the present approach. 

In Sec. 2 of the present paper, after a review of the basic properties of Bogoliubov 
coefficients, we begin the above process of relating them, in our problem, to the complex 
coefficients z n £ of Eq.(3.3) of [3]. Sec. 3 treats particle emission rates. In Sec.4, we consider 
the resulting probability distribution, in terms of the steady-state Bogoliubov coefficients. 
Sec. 5 contains a brief Conclusion. The Appendix relates this probability distribution to 
the Wigner quasi-probability distribution function for harmonic oscillators. 

2. Bogoliubov transformations 

Consider first the ( nearly- )Lorentzian collapse problem. We assume the existence of 
a spherically-symmetrical Lorentzian-signature 'reference' or 'background' metric in the 
form 

ds 2 = - e b ^dt 2 + e a ^ r) dr 2 + r 2 (d6 2 + sin 2 6 d(j) 2 ) . (2.1) 
In this context, it is conventional to define the 'mass function' m(t, r) by 

exp(-a(t,r)) = i _ ^) . (2 . 2) 

In the boundary-value problem outlined above, we write (7^ , <fr) for the 'background' 
spherically-symmetric metric and scalar field, and V M for the background covariant deriva- 
tive. 

Now consider a classical linearised solution <f>^ of the massless scalar wave equation 

V M V M (1) = . (2.3) 

Suppose first that one is given two spherically-symmetric Cauchy surfaces in the space- 
time [14], an initial Cauchy surface S~ and a final Cauchy surface S + . Then, subject to 
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regularity and spatial fall-off conditions, one can expand (p^ in terms either of a basis 
{fu>'£m(x)} of mode solutions Eq.(2.12,16) of [3] adapted to S~ , or of a basis {pu,£ m (x)} 
adapted to S + : 



<a (1) w = E E 

£=0 m=- 

oo £ 

= E E 



=0 m= 



do;' 



dw 



Coj'lm fu'£m(x) + C.C.j , 
&w£m Pu£m(x) + CC. 



(2.4) 

(2.5) 



The c- number complex coefficients {c w /.£ m } are a set of position-independent 'Fourier ampli- 
tudes' labelling the configuration of the field on S~ , while the {b w £ m } refer correspondingly 
to S + . On S~ , the {f w 'em} are an orthonormal, complete, family of complex solutions of 
the wave equation, which contain only positive frequencies (u/ > 0). Here, in addition to 
making use of the spherical symmetry of the background 4-geometry, one adopts a simple 
definition of 'positive-frequency' in regions where the background gravitational and scalar 
fields are approximately static. Similarly, on iS + , the {£W m } are an analogous family of 
positive- frequency solutions (a; > 0). 

In the presence of a future event horizon, H + , the final Cauchy surface S + (say) 
will typically need to cross H + . This is the case in the familiar treatments of black-hole 
evaporation [9-12], which are based on Bogoliubov transformations. Of course, given our 
analytic-continuation strategy, as developed in earlier papers [2,3], with the time-interval at 
infinity taken to be T = \T\ exp(— id) for < 9 < n/2 , one has effectively no event horizon 
in the nearly-Lorentzian problem. Rather, one expects the boundary-value problem to be 
strongly elliptic, and so to resemble the real Riemannian boundary-value problem, with 
no obstruction across any 'horizon'. Thus, in our case, one can replace the usual Cauchy 
surface S + , which in the centre of the space-time lies close to the Lorentzian curvature 
singularity, with the final surface E^ above, on which the final data (hij<p)F are posed. 
One can take the intrinsic geometry of E_p to be nearly flat, with E^ diffeomorphic to M 3 . 
The time-separation from the initial surface E j to E^ , as measured at spatial infinity, is 
T = \T\ exp(— i6), with \T\ very large. A possible choice for a hyper-surface which crosses 
7i + smoothly is a constant 'time' slice in a Painleve-Gullstrand-like coordinate system, 
modified slightly to account for the slow change of the black-hole mass due to emission; 
these coordinates are stationary and non-singular across H + and have been used in a 
tunnelling interpretation of black-hole evaporation [15,16 ]. 

Thus, loosely speaking, we shall consider a 'Bogoliubov transformation to a smooth 
surface, long after the Lorentzian singularity'. In the same context, one may say that 
'the singularity is simply by-passed in the analytic continuation'. That is, there should 
be an analytic complexified classical solution to the boundary-value problem, which only 
reaches a singular boundary precisely at Lorentzian signature (9 = 0). One thus expects to 
have the possibility of circumventing Lorentzian singularities by deforming time-intervals 
suitably into the complex, much as one avoids singularities of functions f(z) in the ordinary 
complex 2-plane by deforming contours. 
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Suppose that the surfaces Ej and Hp are taken to be spherically symmetric. Then, 
the bases above may be chosen to have the form 



<P.im(x) = N(u)^^-e-^YU^) (2-6) 

in terms of the spherical harmonics Yg m (Q) [17], where N(u) is a normalisation factor, and 
where positive frequency corresponds to u> > . On Ej , one chooses the large-r behaviour 
of Ru,e(r) such that w ^ m (a;) is proportional to e~ tuJV at large r , where v = t+r* . Thus, the 
{fuj'em} are purely ingoing at infinity. At E^ , one requires that <j> U £m(x) be proportional 
to e~ tuJU at large r , where u = t — r* ; thus, the {p LU e.m} are purely outgoing at infinity. 
We repeat the reasonable assumption for the collapse problem, that the initial spatial 3- 
geometry and scalar field (the 'star') can be taken to be approximately static. Typically, 
the initial Cauchy surface Ej would be taken to be in the early nearly-static region, whereas 
E p would be at a time long after the collapse, where the background geometry would be 
approximately that of a Vaidya solution with slowly- varying mass [13]. 

The normalisation factor N(u) is determined through use of the natural inner product 

[11] 

(</Wm, <j>uj'£'m') = - i da^ ((pcjem V^* : Vm i - 4>*u'l'm' V^<?Wm) • (2-7) 

For a pair of classical solutions </w m , fyu't'm' ■> this inner product is independent of the 
particular choice of (asymptotically-flat) space-like hypersurface E . The surface element 
da^ in Eq.(2.7) is evaluated in the background space-time. The inner product ( , ) has the 
properties 

{4>u>Zm > ^4>u't'm') = A* (4>ulm > 4>oj't'm') i (2-8) 
(A0o,£ m , 4>u't'm') = A (4>w£m , 4>u}'t'm') , (2.9) 
{4>oj£m , <j>u'e'm')* = {4>uj't'm' , 4>ojim) , (2-10) 

where A is a complex number. For positive-frequency solutions of the massless scalar wave 
equation, this inner product is positive-definite [11]. By a suitable change of basis on each 
of Ej and Tip independently, one can normalise the {fujim} and the {pu,£ m } such that 

tfuim, fu'l'm') = {Puolm , Pu'l'm') = ^££' 5 mm ' S(oj , Oj') , (2-11) 
(futmi fu'l'm') = (Puim 5 Pu't'm') = ® > (2-12) 
(fulmi fu'l'm') = (Pw£m» Pti'i'm') = ~ S W S mm> K u » <*>') . (2.13) 

Since the {f w '£ m } form a complete orthonormal set on Ej , one may expand out a 
typical basis function (solution) p u t m on the surface £/ in the form 



Pujlm 



J <ko' {a^i^im fu'lm + Pu>'u>£m fu'£-m) > (2-14) 
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where the sets of complex numbers {a u '^ m } and {f3 UJ 'ujf.m} give the Bogoliubov coefficients 
[11]. Here, from the definition of spherical harmonics Ye m (Q) in [17], one has Y^ m = 
(— l) m Y^_ m . Of course, the (approximate) spherical symmetry of the background implies 
that solutions with the same I and \m\ are connected in Eq.(2.14). Further, spherical 
symmetry implies that the Bogoliubov coefficients are independent of m , and we shall 
henceforth ignore this index. In the gravitational-collapse problem, the time-dependence 
of the collapse geometry gives f^gm 7^ p u em ; thus, the coefficients f3 LU ' LU g are non-zero and 
there is mixing between positive-and negative-frequency solutions. 

The Bogoliubov coefficients may be expressed in terms of inner products as 

auj>ut = (fWm , fu'tm ) : ( 2 -15) 
Pu'ut = ~ (Puim , fu%-m) • ( 2 - 16 ) 

Given these relations, one may invert Eq.(2.14) by expanding out fui'im in terms of the 
{Putm} and {pZi- m } basis on E F , as 

/•CO 

Jo 

One can similarly relate the 'Fourier amplitudes' {£wm} and {c^'em} above: 

poo 

buiim = / du' {a* UJLoi Cu'im - Pu'ui c u'£,-m) > ( 2 -18) 
JO 

Cto'im = I did [oLu'ul Kim + Ki'ul Kjl-m) • ( 2 -19) 

Jo 

Further, on substituting Eqs.(2.14,17) into Eqs.(2.11,13), one finds that the Bogoliubov 
coefficients must obey the quadratic relations 



I 
I 
I 





00 





00 



du (a u » u t a*, w j - p u , u t PZ,"ut) = u ") > ( 2 - 20 ) 

( a u>ive f3u"uj£ - Puj'ui a u"ue) = > ( 2 - 2 l) 
du" (a U "ui Puj"u>'e - Poj"u>e = . (2.22) 



These equations (2.20-22), which naturally express the property that the Bogoliubov co- 
efficients are the matrix components for a map relating orthonormal bases on Ej and E_p , 
are discussed further in [11]. 

This formalism, involving the use of Bogoliubov coefficients as in the treatment of 
black-hole radiance and other consequences of quantum field theory in curved space-time 
[9-12], can now be connected with the alternative formalism set out in Sees. 2-4 of [3]. 
Consider again the (nearly-) Lorentzian-signature case, where T denotes the proper-time 
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interval between the initial and final hyper surfaces. Making the mode decomposition with 
respect to 'Lorentzian coordinates' (t,r,0,<j)): 



^ oo t 

0«(t,r,0,0) = - E Y £m (n)R £m (t,r) , (2.23) 

£=0 m=-£ 

one arrives at the (£ , m) mode equation: 

(>- a )/ 2 <9 r ) R em - [d t ) 2 R im - -(d t (a-b))(dtRi m ) - V t {t,r)R tm = . (2.24) 
Here, 

m , T) = £^(^ +1) + M^)) , (2 . 25) 

where m(t, r) is defined in Eq.(2.2). For adiabatic modes (that is, for frequencies k of oscil- 
lation which are rapid compared to the time rate of change of the background geometry), 
one has approximately-separable solutions of the mode equation (2.24), of the form 

Rem(t, r) ~ exp(zfct) 6^ m (t,r) , (2.26) 

where £k£mit, r) varies 'slowly' with respect to t . The definition of £fc&n(t, r) is tied down 
completely at spatial infinity (r — > oo), where i^ m (t, r) is required to reduce to a flat space- 
time separated solution, for which £,kem(t,r) loses its t-dependence. By completeness, 
the general solution 4> \t, r , 0, <fi) °f the linear wave equation (2.3) may be written as in 
Eq.(2.18) of [3]: 

= - E E f ^owmW^l^irUfi) • (2-27) 
r to m =-' J -°° sm^jj 



where the {a/^ m } are real coefficients. Here, for simplicity, we are taking, as in [2,3], 
boundary conditions on <f>^ for which = but (j>^\-^ ^ 0. Physically, if the 

gravitational initial data were also taken to be spherically symmetric, this would 

correspond to scalar particle production resulting from isotropic collapse data. At very 
late times, near the final hypersurface E p , the geometry also varies extremely slowly 
with respect to time, and the 'adiabatic separation functions' £,kim{t,r) again reduce to 
functions of r only, namely £fc^ m (r). These obey the mode equation [1-3] 

e^ 2 £ (>-«V 2 ^ + ik 2 - Vt) e« = , (2.28) 

with boundary conditions described below. The spherical symmetry of the background 
implies that £kemi r ) is independent of the quantum number m. Thus, 

iumir) = iteir) . (2.29) 
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Evaluating Eq.(2.27) at the final boundary S F , one finds that 

0(1) ( X ) = -J2 Yl / *^m6fW^m(n) • (2.30) 

Ef r fcO m=-l J -°° 

This relation can be inverted (see Sec. 3 of [3]), given the normalisation of the {£kl( r )}, to 
find the coefficients ak£ m in terms of the final data 4>^ l \x). Indeed, the a^i m characterise 
the final scalar data. Since the mode functions {£,ke(r)} on Ep are real, one can rewrite 
this in the form 

involving only positive /c-values. Here, the geometry in this space-time region is expected 
to be approximated very accurately by a Vaidya metric [13], corresponding to a smoothed- 
out luminosity in the radiation, which varies only slowly with time. Such a metric can be 
put in the diagonal form (2.1), with [13]: 

e~ a = 1 - 2m(t > r) ; e b = ( J^X e~ a . (2.32) 



f(m) 

Here, m(t,r) is a slowly- varying function, with m = (dm/dt), and the function f(m) 
depends on the details of the radiation. 

The 'left' boundary condition (3.1) of [3] reads 

Zki = (2.33) 
at r = . The 'right' boundary condition (3.3) of [3] requires that: 

£ k i(r) ~ (zkt exp(z£x*) + z* u exp(-z/cr*)j , (2.34) 

as r — > oo . A generalisation r* of the standard Regge- Wheeler coordinate r* for the 
Schwarzschild geometry [18] may be defined, by 

° e^/ 2 ^- . (2.35) 



dr* dr 

Here, for each (k,£), z^t is a dimensionless complex coefficient. In our space-time, for 
sufficiently large r , a coordinate r* may also be defined by 

r* = r + 2Mln((r/2M)-lj , (2.36) 

where M is the total ADM (Arnowitt-Deser-Misner) mass [18,19]. For very large r, r* 
and r* are asymptotically equal. 



In the nearly-Lorentzian case, the late-time behaviour of quantum amplitudes will (as 
in Sections 3,4 of [3]) be dominated by the eigen-frequencies k n defined in Eq.(4.1) of [3]: 

kn = y- , (2.37) 

(n = 1, 2, 3, . . . ). Accordingly, let us discretise the frequency integral in Eq.(2.31) above. 
Define real functions {fne(r)}, for n = 1, 2, . . . ; £ = 0, 1, 2, ... , such that 

U{r) = (z nl e ik « r * + z* ne e ~ ik ^) f ne (r) , (2.38) 

where f n £(r) — > 1 as r — > oo , f n e(r) ^ as r — > , and r* ~ r* as r — > oo . 

The discrete description of Eq.(2.31) above is appropriate for the final data (j)^ (x) | E , 
but we further need a discrete description of the space-time behaviour of (f)^(x). For 
this, define Ak n = n/T , and then (at least in a neighbourhood of Hp) define complex 
coefficients A n g m , complex numbers z n e and further real functions {g n e(r)} such that 

= J E Ak ^Anim [Bnee-^-^+Zrue^-^grui^Y^) , (2.39) 

where the {g n i(r)} obey the same boundary conditions as the {f n e(r)} • Eq.(2.39) agrees 
with the discretised version of Eq.(2.31), provided that 

9ni(r) = fni(r) , z n£ = z ne e " ifc " |T| , 

^-n^m = Q>n£m d—n£m 

In order to make a comparison with the Bogoliubov description, consider now the 
discretised version of Eq.(2.5), appropriate to the final surface Hp , taking 

Pntm (x) = N(k n ) e -**»(*-r-) y £m(0) . (2.41) 

Given the boundary condition of regularity at r = , Eq.(2.7) shows that the normalisation 
factor is 

N(k n ) = (27rk n )- 1 / 2 . (2.42) 
Then, comparing the discretised form of Eq.(2.31) with Eq.(2.39), we find 

Pnt(r) = fnt(r) , (2.43) 

z n£ m At\ 

Kim = mi-) • (2 - 44) 

In calculating particle emission rates in the following Section 3, we shall need 

2 2 2 

\b n £m\ = 2nk n \z n £\ \a n £ m + a- n £m\ = 2ft m (k n ) , (2.45) 

9 



where we have used the fall-off property as \k\ — > oo : 

\flm(k)\ ~ \k\~ 3 (2.46) 

of Eq.(4.6) of [3]. For later use, we also record an expression for |£w m | 2 in terms of 
Bogoliubov coefficients, where Eq.(2.18) has been used. Including an implicit summation 
over the index m , we find 



\buim\ 2 = J du' J du" ^a*,^ + PZ"u£ Pu'uiJ c^'im ct"£m 

/•oo />oo / 

-2 J du' J du" Re( P^'^e a*,^ c^'im Cu"i-m 



(2.47) 



3. Particle emission rates 



As remarked in Sec. 5 of [2], the total energy of the final Einstein/scalar field con- 
figuration must equal that of the original configuration, namely M , the true ADM mass 
of the 'space-time'. Recall that we are here considering gravitational and massless-scalar 
perturbations about a given spherically-symmetric background solution, of the form = 

1»v + e h$ + e 2 h$ + . . . , = $ + e <j>W + e 2 </> (2) + ... . The 'second-variation' 0(e 2 ) 
contribution to M (that is, the leading contribution beyond the background contribution) 
on is 

= [ d 3 x(-^H , (3.1) 

where 7i = e~ T^ t and where T^J is the perturbed energy-momentum tensor, given 
in Eq.(3.28) of [2], which is of quadratic order in the scalar-field perturbations, and, if 
evaluated at a late time, has no contribution from the background scalar field <!>(£, r), in 
view of the final boundary conditions that the gravitational data are spherically symmetric 
on Ep. From the wave equation (2.3), integrating by parts, using the boundary conditions 
of regularity and the asymptotic behaviour r 2 <f>^ d r <j)^ = 0(r _1 ) as r — > oo , we find 



£(2) = IJdQj drr 2 e^' 2 </>« 2 - 0« ^ - Q (a - b) 0« 0« 



■ (3-2) 



We may neglect the terms in Eq.(3.2) with integrand proportional to (a — b), compared 
with the [<^>« 2 - 0( 1 >0( 1 )] terms, provided that the adiabatic approximation of Sec. 2 of 
[3] holds; that is, that typical perturbative frequencies u obey u — b\ . 

In the adiabatic case, E^> decomposes into a 'sum' over frequencies. On substituting 
the representation (2.5) into Eq.(3.2) and using Eq.(2.41), we obtain 

/•OO 

This can also be re-written in the 'harmonic-oscillator coordinates' {c^im}, on using 
Eq.(2.47). 



10 



One can read off the final particle-number spectrum from the total energy by writing 



= I duiu^^ . (3.4) 

im 

On using Eq.(2.45), one obtains 



\bcoem\ = 2f £m (uj) , (3.5) 



where fz m (w) is defined through Eqs. (2.41-43). Let AG^ be the number of states or phase- 
space cells in mode u> , and let N' U £ m be the number of (scalar) particles in these states. 
Then 

A/7 

Ari = < n uim > (3.6) 

is the mean number of particles in each of the (quantum) states in mode (oj£m). Since the 
number of states emitted in the frequency interval (u , u> + du>) is 

AC = — du , (3.7) 

TV 

the number of particles in the range (u , u + du>) is 

dN'uim = — < n U £ m > du . (3.8) 

TV 

In the continuum limit, 

dNLtm = < — > doj (3.9) 

TV 

is the total emission rate in the range (u> , u> + cLu), where N^im * s the number of particles 
emitted over the duration \T\ . (The factor tv , rather than the usual factor 2tv , occurs in 
the denominator in Eq.(3.9) due to the initial time being t = rather than t = —\T\ .) 
Alternatively, Eq.(3.9) just gives the number of particles passing out through the surface 
of a sphere centred on the collapsing 'star' or black hole. 

In the case (as here) that the particles emitted do not decay further, one identifies 
< n u tm > as [11] 

<n U £ m > = \M , (3.10) 

where \(3 u i\ is defined by 



du;' fc ue fc u »i = \PJ 2 S(u,uj") . (3.11) 



POO 

Jo 



'0 

Then the total number spectrum, added up over the black-hole life-time, is 
dN'„ 1 ^ „ f tF 



I ^ 



E ' -dM 



dt 



(3.12) 
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where N S0J e m counts the number of states with spin s and quantum numbers (u>£m) . In the 
original calculation of Bogoliubov coefficients and of probabilities for particle emission by 
a non-rotating black hole [10], it was found (neglecting the effects of back-reaction) that 

l/^fcnl 2 = T autm (m) (e 4 ™-(-l) 2s ) -1 , (3.13) 

where T 8UJ £ m (rh) is the transmission probability over the centrifugal barrier [20] and m = 
2Muj is dimensionless. This calculation, of course, referred to the case where one did not 
have a final surface of topology M 3 . But, because of the very-high-frequency (adiabatic) 
method through which this expression for \(3 8 uim\ was calculated, it should still be valid 
(up to tiny corrections) in our present M 3 case. 

The semi-classical expression for the rate of mass loss is [21] 

dM q (M) 

dt M 2 ' 1 J 

where cko(M) ~ constant for massless or ultra-relativistic particles, and M = M(Mj,t) is 
the mass to which a black hole of initial mass Mj has been reduced after time t . Turning 
again to the total number spectrum over the black-hole life-time, given by N' ' w as in 
Eq.(3.12), note that, for massless particles, the Bogoliubov coefficients \(3 SLJ £m\ depend on 
u> only through m (on dimensional grounds). Hence, 



dN' u 1 



_ /•2Mjai 

/ dm m 2 \(3 SOJ t m (m)\ 2 . (3.15) 

_,.„ Jo 



duo 8ncto u 3 

slm 

In particular, in considering the high-frequency limit ujMj > 1 of this expression, we 
may replace the upper limit in the m-integral by infinity. From this, one finds that the 
time-integrated energy distribution of all the massless particles radiated by the black hole 
(of initial mass Mi) is 

— ~ Cl u , (3.16) 

for uj ^> (M/) _1 , where c\ is a real number, so giving the form of the high-energy tail 
of the spectrum. However, most of the particles produced during the quasi-stationary 
regime of the radiating black hole, before its expected explosive phase with M — > , are 

at a temperature corresponding to the initial mass of the hole (uj ~ (M/) _1 j, where the 

energy distribution peaks. 

How is the above analysis related to our theory which includes back-reaction? One 
might expect that the high-energy behaviour of the {a w £ m } 'coordinates', which describe 
the perturbed scalar field (p^ on the final surface Tip , would be related to the emission in 
the final moments of evaporation. This can be seen from Eq.(3.5), taking the frequencies 
k n = &n = mr/\T\ , without an explicit calculation of the Bogoliubov coefficients. The 
high-energy behaviour of the function fe m (uj), which is defined through Eq.(4.3) of [3], was 
found in Eq.(4.6) of [3] from the requirement that the contour at infinity in the integral (4.2) 
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of [3] for the action should contribute zero. For massless scalar particles, this requirement, 
expressed through Eq.(3.5), reads 



oc K) , (3.17) 
for a n ^> (M/) _1 , in agreement with the high-energy behaviour in Eq.(3.16). Here, 

(2) 

again, the requirement for a finite imaginary part for the classical Lorentzian action S^ ss , 
corresponding to a non-trivial probability distribution over final configurations, is linked 
with the behaviour of the high-energy number spectrum; the latter is of course subject in 
principle to observational test. 

It is well known that the thermal equilibrium between a black hole and the exterior 
radiation is unstable, as the specific heat in the canonical ensemble is negative [22]. The 
canonical ensemble breaks down for black holes since the canonical partition function di- 
verges for all temperatures. Rather than work at fixed temperature, one must consider 
the micro-canonical ensemble, which is tailored to configurations of fixed energy. At the 
high-energy end of the emission spectrum, therefore, when the black hole approaches the 
Planck scale, the canonical distribution, as given by Eq.(3.13), must be replaced by the 
micro-canonical distribution. Naturally, the mass-loss rates for the two ensembles differ 
considerably. The micro-canonical decay rate modifies the small-mass behaviour of the 
black hole, where Eq.(3.14) breaks down. For the low- frequency quanta {uj <C M) char- 
acteristic of the majority of the evaporation process, the canonical and micro-canonical 
ensembles are almost equivalent, and one obtains a Planck-like number spectrum Eq.(3.13) 
[23]. 

Prior to the black hole's total evaporation, it is possible for the black hole to emit a 
single high-energy quantum with energy comparable to the initial black-hole mass. As \(3j \ 
is the average number of quanta present, then in the high-energy tail of the spectrum, one 
has \j3j\ ~ e , where S is the Boltzmann entropy of the system, since there is only one 
state out of a total of e states for which all the energy is concentrated into one quantum 
[24]. More precisely, 

1/3 I 2 

~ ASbh , (3.18) 



1/3/ 



l«jl 



where ASbh is the difference in black-hole entropy before and after emission. This result 
is expected to be valid for a general spherically-symmetric black hole. 

4. Probabilistic interpretation 

When \T\ is large but finite, one finds from Eq.(4.16) of [3], that the density function 

P[{<W;|T|] = Ne- 5 ^ M < exp(-2Im 5^[{a fc/m }; |T|]) , (4.1) 

where N is a suitable normalisation factor, describes the conditional probability density 
over the final perturbative scalar boundary data , the condition being that the 

perturbations obey the initial conditions = of Eq.(2.7) of [3]. The limit 5^0 
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should then be taken; the only data at spatial infinity itself consist of the proper time |T| . 
This probability is, of course, the squared norm of a complex quantum amplitude, in the 
context of [3]. 

Even though the probability distribution (4.1) arises from squaring a quantum ampli- 
tude, one can still ask, given the thermal nature of black-hole evaporation in the usual de- 
scription, whether the probability distribution (4.1) can be viewed in terms of the diagonal 
components of some non-trivial 'density matrix' in a suitable basis [9-12,32]. To evaluate 
more explicitly the probability distribution, consider Eq.(2.36), using also Eqs. (2.20,45,47) 
and Eq.(4.16) of [3]. One finds, in the continuum limit of large |T|, taking the {aum} as 
the final 'coordinate' variables, that: 



ujlm 



N expl -2V / du \b, 
iVexp^-2V / du \cutm 



Im 
"OO 



+ 



poo poo poo 

2 du' du" Cu'tm cim m \ du Pu"u>£ Pu'ui 

Jo Jo Jo 



du' 



du Re I c^'Cm Cojiig^-m 



du (3 i0 n Lol a* w 



(4-2) 

The first two terms inside the square bracket are positive-definite, whereas the third is of 
indefinite sign. The situation becomes much clearer in the case that 



duJ Pu"u>£ & 



l/^UI 2 5{u\u") , 



/•OO 

/ du at'ue = . 

Jo 



(4.3) 



(4.4) 



This holds, in particular, for the steady-state Bogoliubov coefficients in the calculation 
which neglects back-reaction on the metric [10]. It remains to check whether this diago- 
nal form persists under adiabatic propagation through a slowly-varying potential. From 
Eq.(2.20), one has 



\<Xu'i\ 



u>'£\ 



= 1 



(4.5) 



Thence, 



P 



= N exp 



/•o. 



du \c 



wim | 



£m 



du \/3 ui \ | 



(4.6) 



N ] [ exp 



nlm 



2 2 

-2 (Au n ) \c nim \ - 4 (Au n ) \Pni\ \c n£m \ 
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where Au n = ir/\T\ . The product over ntm tells us that the modes evolve independently; 
such uncorrelated modes are a consequence of the linearised theory. 

A further argument for the diagonal form of the Bogoliubov coefficients is through 
the tunnelling interpretation of black-hole evaporation. Particles are created in pairs just 
outside the future horizon 7i + , with one member always falling into the singularity whilst 
the other escapes to infinity, or is reflected back down the hole. By constructing a ba- 
sis of scalar field modes which is continuous across the future horizon, we automatically 
incorporate both the positive- and negative-energy particles [25]. To achieve this in the 
Schwarzschild picture, however, one must avoid the coordinate singularity at the future 
horizon by adding a small imaginary part to the mass M. Such a procedure is in keeping 
with Feynman's complex-time technique - which we expound in this paper - since the con- 
jugacy between the intial mass Mj and asymptotic time T implies that moving T slightly 
into the lower complex plane is equivalent to adding a small positive imaginary part to 
Mj. The Bogoliubov coefficients are just the weights of the outgoing and ingoing particle 
components of the scalar field across the future horizon. 

Equation (4.6) expresses the probability distribution of final field configurations in 
terms of amplitudes associated with the S~ Cauchy surface. The nature of Eq.(4.6) sug- 
gests that there is a shift in the width of the ground-state probability distribution. That 
is, -P[{c n ^ m }] appears as an excited state relative to the ground state probability density 
P[{b n £m}]- Having averaged over {cj} amplitudes, the stability of the S~ state is guar- 

2 2 

anteed by the fact that \aj\ + \(3j\ > for all j. However, the shift only occurs for 
the infra-red or low-frequency components of the field, if we assume that the Bogoliubov 
coefficients \(3j\ 2 decay (exponentially) rapidly at high frequency and are finite for low fre- 
quency. This is just an expression of local (coordinate) covariance, in that the T,p and S~ 
representations are equivalent for the ultra-violet properties of the theory. An infra-red 
ambiguity is manifest through the occurrence of the Hawking effect with non-zero \/3j\ . 
Hence, the Hawking effect is a consequence of the ambiguity in the choice of space-like 
evolution hypersurfaces [26] . 

5. Conclusion 

In this paper, we have made a connection between the present boundary-value ap- 
proach, as used in this study of quantum amplitudes in black-hole evaporation, and the 
original approach by means of Bogoliubov coefficients. This connection is established 
through consideration of the radial equation obeyed by the (adiabatic) perturbations; in 
this paper, it is the spin-0 perturbations which are studied. The discussion in Sec. 4 gives 
the probability distribution for configurations of the perturbative scalar field, on a final hy- 
persurface Sj? at a very late time, by taking the squared norms of amplitudes, as a product 
of independent Gaussians. Since our approach gives a final pure state, whereas the tradi- 
tional approach gives a non-trivial density matrix, one might expect some difference (not 
necessarily great) between the probabilistic predictions of our boundary-value approach 
and of the Bogoliubov-coefficient approach, in the case of gravitational collapse to a black 
hole. In the Appendix, we consider a somewhat different question, and find a possible 
interpretation of our probability distribution, in terms of the Wigner quasi-probability 
distribution for harmonic oscillators. 
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Further work in this area has involved the study of approximately Vaidya-like metrics 
[13], as giving an accurate approximation to the gravitational field in the region of space- 
time containing the flux of outgoing black-hole radiation. At the same time, we discuss 
the radiative spin-0 (scalar) and spin-2 (graviton) fields which both act as a source for 
the gravitational field and propagate within it. This Vaidya-like description is essential 
in the treatment above of the adiabatic perturbation modes. In later work, we shall 
give a further alternative description of the quantum states found in our boundary-value 
approach, namely, a description in terms of coherent and squeezed states. In [27,28] we 
repeat the spin-0 boundary- value quantum calculation of [3] (as used in the present paper), 
but for other spins, including the fermionic spin-| (neutrino) case. 

Appendix: Density-Matrix Interpretation and Wigner Distribution Functions 

Our probability distribution -P[{GWm}] = -P[{c^ m }] above naturally arises from the 
pure state that we have discussed. But we now check whether or not -P[{cw m }] can 
be expressed in terms of the diagonal elements of some other kind of 'density-matrix 
distribution'. One might expect that a probability density for the final scalar field would 
incorporate not only the randomness in the particle-emission process but also possible 
choices for the boundary data. In light of the loss of phase information in Eq.(4.6) with 
respect to the {c^im} amplitudes, a diagonal density-matrix-like expression for -P[{cw m }] 
might now be possible. Eq.(4.6) can be rewritten so as to make this clear, when one 
employs the generating function for the Laguerre polynomials {L k (x)} [29]: 



1 



(1-8) 



exp 



— xs 



= L k( x ) sk > N < 1 



k=0 



where 



Here, we set 



L k(x) = Yl 



n=0 



k\ (-x) 7 



s = 



(3, 



nl 



OL n l 



nl nl 



< 1 



and 

x = 2X n £ m = 4 (Au n ) \c n £m\ , 
a dimensionless quantity. Then, from Eqs.(4.6,Al), and writing j = nim, we have 

{ Cj }] = N]J expi-Xj) exp (-2 A,) 

3 

oo 

= N J] exp(-A J ) P ^j) L kA 2X j) 

3 kj=0 

= N P ( k 3) II ex P("^) LksPXj) , 

{*,•} 3 



(Al) 

(A2) 

(A3) 
(Al) 



(Ah) 
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where P(kj) is denned to be 



1 



a. 



Pi 



a 3 



P(kj) = --2 ^ , (^6) 



namely the probability to observe the field in the state fc, such that from Eq.(2.20) one 
has [30] 

oo 
fc a -=0 

Alternatively, P(kj) is the probability of finding particles outgoing at X + (future null 
infinity), in the mode labelled by nlm . The average number of particles outgoing at 1+ 
per unit frequency around u n per unit time is then < kj >= \(3j\ , independently of 

m . The factor ^ 2 effectively controls the probability per unit time to emit particles in 
the frequency range around ui n . Each radiation state is uniquely specified by the set of 
occupation numbers {kj} of the various modes {n£m} and occurs with probability P({kj}). 

Eq.(A5) describes the probability distribution -P[{cj}] for the final scalar field, with 
the help of the Bose-Einstein distribution (A6), which describes the randomness in the 
particle emission process. Since, from Eq.(4.6), P[{cj}] depends only on the modulus 
\cj\ of the complex number Cj , the sum in Eq.(A5) is diagonal. Of course, the explicit 
expression for the 'probability' (A5) for a given total mass Mj is clearly complicated, since, 
as in Eq.(3.13), the Bogoliubov coefficients \fij\ depend on the transmission probability 

ry 

Before we investigate further the probability distribution (equivalently -P[{cj}]),| 

we define, for fixed j , the measure d 2 cij on the complex aj-space to be the natural measure 
on R 2 : 

d 2 a = d[Re(aj)] d[lm( aj )] . (A8) 

Since the probability density only depends on \a,j \ , it is further convenient to rewrite the 
measure in terms of polar coordinates, through a,j = Rj e %di , for each j . We then choose 
the normalisation factor N in Eqs.(4.1,2,6,A5) such that 



\ f II d 2 ajP[{aj}] = 1 . (A9) 



3 



In the language of Eq.(4.6), this involves a Gaussian integral over the infinitely many 
variables a,j or Cj , which will require regularisation, for example by the (- function technique 
[31]. This procedure will be simpler when one works with a locally-supersymmetric theory, 
such as the theory of N = 1 supergravity with massless-scalar/spin-^ supermatter [2]. 

In that case, there should not be any divergences in the calculation of N, because of 
cancellations between the bosonic and fermionic degrees of freedom, and such cancellations 
should make the detailed expression for N much simpler. 

In Eq.(A5), the Laguerre polynomial Lfc(x) has k real, distinct roots; further, for 
positive argument Xj > and excited states (kj > 0), the function e~ Xj L^ j (2Xj) takes 
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negative values for certain ranges of Xj . This function cannot therefore be interpreted 
as a probability density. Regarding the question of whether the probability distribution 
P[{aj}] might arise from the diagonal components of some kind of 'density matrix', let us 
consider an Hermitian position-space density matrix of the form 

p(x,y) = w * ^0*0 V£(y) , (^10) 

i 

where the {ipi(x)} form a complete set and Wi is the probability to be in the state ipi(x). If 
the system is in the single state defined by ifii (x) , then the probability density for observing 
x is . The full p(x, y) does not have a probabilistic interpretation, but the diagonal 

contribution 

P{x) = p{x,x) = w * \M%)\ 2 (All) 

i 

is the probability density for observing the coordinate x . For a pure state, one has 

p(x,y) = . (A12) 

Therefore, P(x) for a pure state is the squared norm of a complex quantum amplitude. 

In our case, if the field is in the k-th state, then according to Eq.(All), N e~ Xj (2X,)| 
would be the 'probability density' for the 'harmonic oscillator coordinate' Xj . Although 
this function takes on negative values, the Bose-Einstein distribution Eq.(A6) effectively 
smoothes out the oscillations in Ne~ Xj (2Xj) to give a positive probability density 
when we sum over all possible states {kj}. Although this provides an interesting way of 
breaking down the probability distribution (A5), it does not provide us with a sensible den- 
sity matrix, derived post hoc from the quantum amplitude exp^z J S , ^ ss [{a^ m }; T]j , where 

S SL is S iven b y Eq.(AlO) of [3]. 

The density Ne~ Xj L k . (2Xj), however, is closely related to the Wigner quasi-probability| 
distribution for the harmonic oscillator [32]. Wigner functions are quantum-mechanical 
analogues of classical phase-space distributions, which are formally equivalent to more fa- 
miliar Hilbert-space or density-matrix formulations of quantum mechanics. However, as 
the uncertainty principle forbids the existence of a probability distribution for simultane- 
ously well-defined conjugate variables, such as coordinates and momenta, Wigner functions 
do not have all the properties of a classical phase-space probability density. In particular, 
Wigner functions must become negative in some domain of phase space. 

For a system in a mixed state described by a density matrix p(x, y) , the one-dimensional^ 
Wigner function is defined as [32] 

From Eq.(AlO), 

W(p,q) = ^T^V^cz), (AU) 
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where 



For a pure state, 

The ^-dependent function W(p, q) is manifestly real and has the properties 

/CO 
dp W(p,q) = \V(q)\\ (A17) 
-00 

/CO 
dq W(p,q) = \V(p)\\ (A18) 
-co 

/CO 
dpdg W(p,g) = 1, (A19) 
-co 

where ^(p) is the Fourier transform of For the fc-th state of the harmonic oscillator 

with the wave function 

1 uj ~ 

i) k {x) = -j= (J^y exp(-ivx 2 /2h) H k (xy/u/h), (A20) 



then 



W k (p,q) = 2(-l) k e W (-2H(p,q)/hu J ) L k (^^y (A21) 

where H(p,q) = (1/2) (p 2 + u 2 q 2 ) is the Hamiltonian and u> is the oscillator frequency. 
Further, using Eq.(Al), W k (p, q) is normalised according to 

CO 

W k (p,q) = 1. (A22) 



k=0 



As H increases, W k (p, q) oscillates with a decreasing amplitude and with increasing broad 
peaks. Introducing the complex variables {cij} which represent a phase-space point (j>j, qj), 
with a mode index j for the infinite-dimensional phase space, where 

aj = J2^ ( Ujl/2qj+lUJ ~ 1/2pj )' {A23) 
one can express W k (p, q) as 

W k ( aj ) = 2(-l) fc e- 2 l^l 2 L fc (4|a J | 2 ). (A24) 

Using Lfe(O) = 1, one has 

W fc (0) = 2(-l) fc , (A25) 
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which is the maximum value of Wk(a,j). Thus, with the replacement 2\aj\ — > Xj , the 
final term in Eq.(A5) can be written as 

1 * ~ W k .(Xj) 

Ne- x iL kj (2Xj) = l NW ki (X j )W kj (0) = N ^ ^ . (A26) 

We can interpret the product - or ratio - of two Wigner functions in the following 
way. Consider two statistically independent subsystems A and B. As in statistical physics, 
the Wigner function W describing the system A U B is given by the product of the Wigner 
functions for the two subsystems 

W(pA,qA,PB,qB) = W A (pA,qA) W B (pB,qB)- (A27) 

Thus, the state of the subsystem B does not affect the 'probabilities' of various states of 
the other subsystem A. Such a factorisation would not be possible if the subsystems were 
statistically dependent, in which case, the Wigner function of each subsystem is defined 
by integrating over the other system's variables: 

W A (pA,qA) = J dq B J dp B W(p A , qA.PB, ?b). (-428) 

In our case, one might interpret the system B, say, as corresponding to the weak initial 
state of the perturbations at early times, and A as corresponding to the final state of the 
fluctuations at late times. Thus, referring to Eq.(All), the iV'iWI term is replaced by 
WkjiX^/WkjiO), which gives a Wigner function for the system A at late times. 
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